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Abstract 

We investigate the effects of chameleon scalar field to the effective density and pres- 
sure of a dark matter halo. The pressure is generated from the chameleonic fifth force 
on the matter. We demonstrate that the thick-shell non-singular boundary condition 
which forbids singular point leads to extremely stringent constraint on the matter- 
chameleon coupling when applied to galaxy. We argue that chameleon profile with cen- 
tral singularity is more likely to develop in general physical situation. The chameleonic 
fifth force from the chameleon profile with central singularity experienced by the dark 
matter could significantly modify the rotation curve of galaxies. The chameleonic fifth 
force could generate steeper cusp to the rotation curves in any dark matter profiles 
starting from the Navarro-Frenk- White (NFW) to the pseudo-isothermal (ISO) profile. 
Upper limits on the coupling constant between the chameleon and the dark matter 
are estimated from observational data of the late-type Low-Surface-Brightness galax- 
ies (LSB). It is in the order of /3 < 10~^. 
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1 Introduction 



With series of observations in the recent years, we come to reahze that the expansion of 
our universe is accelerating [H [21 [3] . However, the exact cause of the acceleration is not 
determined for certain. The popular explanations are based on dark energy models ranging 
from dark energy in the form of the cosmological constant, quintessence, phantom, to the 
scalar-tensor, modified and f{R) gravity models. The chameleon dark energy is one of the 
scalar-tensor theories that use scalar field to drive the accelerated expansion of the universe. 
This model is based on the hypothesis that a scalar field couples with matter via a conformal 
transformation. As a consequence, mass of the scalar field depends on the matter density 
in a significant way. We call this kind of scalar a chameleon scalar field after its ability 
to adapt itself to the environment. Its mass becomes large when the matter is abundant 
and becomes tiny in the low density region. Therefore it can effectively hide itself from any 
detection in such environment [H El El [7]. The effects of the chameleon scalar field on the 
earth are consequently suppressed and the constraints on the fifth force could be evaded up 
to the scale of the solar system. 

Since chameleon can adapt its mass very well to the environment, it can evade most 
gravitational constraints used to constrain other scalar gravity theories. In the models where 
chameleon interact with photons, the most stringent constraints on the chameleon coupling to 
the photon comes from Sunyaev-Zel'dovich (SZ) effect in the Cosmic Microwave Background. 
The SZ measurements of the Coma cluster place strong constraints on the photon-chameleon 
coupling. The bounds are approximately ges < 10~^ — 10~® GeV~^ [8j. Weaker constraints 
come from the Earth bound experiments such as GammeV which excluded the region 2 x 
10"^ < ^eff < 4 X 10-6 GeV-i [gj. 

On the other hand, the strongest constraints on the chameleonic matter coupling comes 
from particle colliders and it is merely /3 < 2 x 10^^ [TUl E] (definition of (3 is given in 
Section 2). This constraint is simply the bound on the new physics mass scale typical for 
any new physics scenario investigated at the particle colliders (Mnew physics > 100 GeV). In 
contrast to the chameleon-photon coupling, no Earth-bound experiments to constrain the 
chameleon-matter coupling is operating. Future tests of gravity in space to measure the 
variation of the gravitational constant could place certain bounds on the matter-chameleon 
coupling dans]. 

A useful insight regarding nature of the chameleon is the fact that even though the 
chameleon field increases its mass to hide its gravitational effect, the chameleon undergoes a 
spatial variation in doing so. When a field varies in space, it generates a density and pressure. 
The gradient of pressure then induces a pressure gradient force subsequently. This pressure 
gradient force is actually the fifth force generated by the chameleon-matter coupling [B] , the 
very interaction responsible for changing the chameleonic mass according to the environment. 
In this article, we will use the effect of the chameleonic fifth force on the matter in the dark 
matter (DM) halo to establish an extremely stringent constraint on the chameleon-matter 
coupling (3. 

The stringent constraint on the matter-chameleon coupling from the chameleonic fifth 
iPor the interaction C ^ -\BF{(j)/M)F^^Ff''' .g^^f f = (lnBF),,/,(0min/M). 
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force is a direct result of the singular chameleon solution inevitably developed in a sufficiently 
large massive object such as the galactic DM halo. Rapidly raising chameleon field near the 
singular point at the center induces appreciably large fifth force to the matter in the galaxy. 
The chameleonic fifth force significantly reduces the circular rotation velocity of matters in 
the core region. It will make the cusp of rotation curves steeper. For the late-type Low- 
Surface-Brightness (LSB) galaxies, the dominant gravitating element is the DM halo. We 
use recent observational data on the rotation curves of the LSB galaxies to place constraints 
on the chameleon-matter coupling. The upper bounds could be as low as /3 < 1 x 10~^ 
depending on which type of DM profile is used. 

This article is organized as the following. In Section 2, we review briefiy on the chameleon 
model and derive the equation of motion governing its profile in the presence of matter. In 
Section 3, we discuss the inevitability of the chameleon profile with singular point at the 
center of the mass distribution. Enforcing the non-singular boundary condition leads to ex- 
tremely stringent constraint on the matter-chameleon coupling /3 (of order 10~^). We argue 
that this boundary condition and consequently this constraint is physically unreasonable. 
The chameleon profile with central singularity is then numerically obtained for general sit- 
uation and the density, pressure and the fifth force in the presence of a chameleon scalar 
field are subsequently calculated. Effects of chameleonic fifth force to the rotation curves 
of LSB galaxies are explored in Section 4 and 5. Constraints on the upper bound of the 
matter-chameleon coupling constant are obtained and discussed in Section 5. Changes of 
the power index in the power-law self potential of chameleon are shown not to affect most 
results. Section 6 concludes our work. 



2 Effective Density, Pressure and the Chameleonic Fifth 
Force 

A chameleon scalar field interacts with matter through a conformal coupling which could be 
absorbed into the matter action as the following 



^R-\{d<pf-V{<p) 



where ^^j, = A'^{(l))g^y, A((f)) = e^'^^^'^^'- is a conformal coupling. 

In order to find the effective density and pressure, we need to determine the profile of 
the chameleon in dark matter (DM) halo. We set the dynamics of the chameleon in the 
thick-shell regime |H HI]. Namely, we assume that the value of scalar field which minimizes 
the effective potential {(pmm) only stays at the exterior of the dark matter halo. The value 
of the scalar field in the interior of the halo is not (pmin but 0(r) which is determined by the 
matter density of the halo. The equation of motion of the chameleon scalar field is [S] 

V V = V,^ +a^PmA (0) , (2) 

where is the chameleon scalar field, pm is the matter density, and = We can use 

available density profiles of dark matter such as the Navarro-Frenk- White (NEW) PJl] and 
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the pseudo-isothermal (ISO) profile [TB] to simulate the effects of the chameleon pressure to 
the DM halo. 

For simplicity, we assume that the dark matter halo has spherical symmetry and the 
chameleon scalar field is static. The effect of spacetime curvature for the DM halo on the 
chameleon profile is negligible. Therefore, Eqn.([2]) is reduced to 

d'^(j) 2d(f) ,,,,,, 

First we will consider the scalar self potential in the form of the inverse-power-law po- 
tential ^(0) = ^^n" jI3 HI HI E]. Even though this form of self potential is extremely 
constrained if not already ruled out by the Lunar Laser Ranging (LLR) experiment and cos- 
mological constraints Ref . [19] , it contains minimal amount of parameters and serves as the 
simplest chameleon model. A more viable self potential of the form = M'^(l-F/i(M/0)") 
which has not been ruled out by the same analyses, having one extra parameter /i, will be 
compatible with our result with the substitution M^+" — )■ /iM^+" (since only V, (p appears in 
the equation of motion). With the power-law self potential, the equation of motion becomes 

The right-hand side of the equation can be defined to be the derivative with respect to 
the chameleon field of an effective potential 

Veff = V + p^{r)e^^""-K (5) 

The value of chameleon (pmin which gives minimum effective potential is then given by 



n 



(6) 



For our purpose, we will define (pmin for Pm = Poo (average density of the universe) ~ 10"^^ 
kg/m^. 

To consider gravitational effects of the DM halo and the chameleon, we solve the Einstein 
equations for a spherically symmetric metric 

ds^ = ~A (r) + B (r) dr'^ + r'^dO'^ + sin^ Odcj)'^. (7) 

The Einstein equations are 

^ + 1^ = (8) 

-8.gt;. (9) 
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The energy-momentum tensor in the Einstein equations are the total energy-momentum 
tensor {T^i^^^i-^ = T^^^itf^y.-) + ^(^p due to the couphng with matter. From action of the 
chameleon scalar field , we obtain 

Tj^i^) = dyd,^ - 6i: (^g'^^d^cpdp^ + V{(P)\ (10) 



and 



T^m = -Pw = - (11) 

(12) 

= P^. (13) 



rpr 








2B ~ 














rp4> 


~2B ~ 





This is an isotropic distribution with curious behaviour which needs extra caution since 
and P^, P*^ are not necessarily equivalent Substitute into the Einstein field equations, the 
metric then has to satisfy 

(P- 1) B' ^( , , 0" M4+"\ , ^ 



Pr2 P2r V 2P 

'^-'^ ^' (16) 



Pr2 rAB \ 2B 

In the equation of motion of the chameleon scalar field, we will set the pressure of matter 
to zero by assuming that the matter in the universe is in the dust form and ignoring the 
possible annihilation pressure of dark matter and such [20] • The effective pressure thus only 
comes from the scalar field. The effective density and effective pressure are then 

Pe// = + ^ + (16) 

iVP^+" 

respectively. Observe that the contribution of /3 comes in the determination of the chameleon 
profile Eqn. (|1]). The effective pressure from the chameleon will thus change with varying /3 
through the changes in the chameleon profile in the dark matter halo. However, the dynamics 
of the matter with respect to the profile of the chameleon can be obtained directly from the 
conservation of the energy-momentum tensor 

V.T(l^,«e.) + V.Tf,'; = 0, (18) 



^In terms of the pressure gradient force when reduced to the Euler equation of the chameleonic fluid, 
there will be extra term ^{2P^ — — P'^) in addition to in the radial direction. This term is required 
in order to obtain the correct fifth force expression as in Eqn. (j20p . 
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which leads to 



or 



pmdtv = --^p„V0, (19) 
IVIpi 



^ V0, (20) 



PI 



by using the equation of motion, Eqn. (jlj), and assuming a pressureless matter (see also the 
Appendix [A]). Namely, the chameleon- matter coupling induces the fifth force acting onto the 
matter (Ref. [6]). This fifth force could change the rotation curve of the galaxy substantially 
provided that the variation of the chameleon is sufficiently large as we shall see later in 
Section 4. 



3 Constraints on the matter-chameleon couphng and 
singular solutions of the chameleon 

In this section, we will demonstrate starting from the thick-shell regime, that the chameleon 
profile within a sufficiently large massive object can satisfy the non-singular boundary and 
positivity condition, 0'(O) = 0, 0(r) > only when (3 < /^max- 

For a region with sufficiently large matter density pm, an excellent approximation of the 
chameleon profile can be obtained analytically. Since Mpi is large (~ 2.4 x 10^^ GeV) and 
M is small (~ 10~^ eV), we can approximate a0p(r)A(0) ~ and neglect the potential 

pi 

term in the equation of motion to obtain 

- (21, 



For a boundary condition 



'"max) 4'mini max 



fmax) — . ,^ 2 ) (22) 



max 



the general solution can be written as 

0'(r) = 4;^WO-^o + 7), (23) 

where 7 represents (proportional to) the slope of chameleon profile at the boundary of the 
mass distribution and Mq is the total mass of the object (e.g. galaxy) at Vmax- This expression 
can be integrated directly for any mass profile M(r) to obtain the corresponding chameleon 
solution. 

There are three classes of positive solutions categorized by the value of 7, 
1. 7 < Mq ; singular at r = 
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2. 7 = Mq ; 0'r^|j.=o = (nonsingular) 

3. 7 > Mq ; truncated at finite r (unphysical) . 



Varying the slope 7 > 0, various chameleon solutions (using the full equation of motion, 
Eqn. dl]), with the potential included) for sufficiently small /3 can be obtained numerically 
clS dbTQ shown in Fig. [TJ All three kinds of solutions are presented. 



U5005 ISO Halo, /3 = 1.69 x 10" 



600000 
500000 




500000 



300 



U5005 NFW Halo, /5 = 1.76 x 10"^ 




Figure 1: The chameleon solutions for various 7. The red curve is the approximated ana- 
lytic solution neglecting the potential V{(j)). The ISO solution (left) satisfies non-singular 
boundary condition 0'(O) = but not the NFW (right). The other curves are obtained 
numerically. 

The third kind of solution is unphysical since it is truncated at finite radial distance. 
Among the remaining physical solutions, the only nonsingular one satisfies the boundary 
condition 0V^|r=o = (class 2). The nonsingular solutions of the chameleon are investigated 
by most literature with the boundary condition (f)'{0) = and finite 0(0). In the following 
subsection, we will demonstrate that for certain matter profiles, this boundary condition 
cannot be satisfied. And for those matter profiles, the most natural chameleon solution 
would be the singular one. 

3.1 Analytic solutions for small potential and remarks on thick- 
shell boundary conditions 

In this subsection, we consider analytic solutions of the chameleon coupling to various DM 
models: the NFW, ISO, and the parametrized model (PM). The potential term V{(f)) will be 
neglected. Subsequently, both the approximated analytic and numerical (with the potential 
included) solutions will be compared and justified the validity of the approximation. We 
will investigate whether which DM profile allows non-singular chameleonic solution with fiat 
boundary condition 0'(O) = 0. 

For nonsingular case in class 2, we have the gradient of the scalar given in terms of the 
matter density as 

0'r2 = j^J p{rydr, (24) 
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where p{r) is the matter distribution. 

Po 

1(1 + r 



For the NFW profile : Pnfw{i^) = tt^-^-^tt^- Substitute into Eqn. ( I2l|) gives 



<P'r' = 4^f^ + ln(a + r))+Ci. (25) 



Mpi \a + r 

Taking hmit r — )■ and assuming 0' — )■ oo slower than in this limit, the constant Ci 
can be determined to be 

Ci = -^(l + ln(a)). (26) 
Mpi 



Therefore 



r = — — — h ln(a + r) — — (1 + ln(a)) 



Mpi \a + r ' V 



^'w = #§(m + i)-^), (27) 



Mpir^ V 0, a + r 

This is the solution of 0' for arbitrary initial condition, 0'(r = 0). The solution of the 
chameleon scalar field in the NFW dark matter halo can then be obtained, 



a'Ppofl ln(l + ^; 



Mpi \a 



+ 0(0). (2J 



Thus, we can see that at the origin of the NFW dark matter halo (r = 0) the chameleon 
profile cannot be fiat with 

.'(0) ^ ^. (29) 

This proves conclusively that for the chameleon coupling to the NFW DM in the thick- 
shell regime, the chameleon field cannot satisfy the boundary condition 0'(O) = 0. Figure [2] 
shows both the approximated analytic and numerical (with potential included) solutions, 
the difference is minimal. 

For the ISO profile: piso{r) = (j_\2 ^ 

/5 



leading to 



(l>'r = ^ / pisoir)r dr, 
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U5005 NFW Halo, /? = 1.76 x lO"'' 
600 000 r 




R(kpc) 



Figure 2: The chameleon profile in the NFW DM halo of U5005 galaxy for /3 = 1.76 x 10"^. 
The analytic and numerical chameleon solutions are represented in red line and black line 
respectively. 



Taking limit r — > and assuming 0' — )■ oo slower than in this limit, we obtain Ci = 0. 
Therefore 

^'M - '''' 

Integrate to obtain solution of the chameleon scalar field in ISO dark matter halo, 

_ /arctan(./fl.) + R-) ^ J_\ 

' Mpi \ r 2R. rJ ' ^ ' 

Thus, the chameleon profile at the origin of the ISO dark matter halo can be flat with 
(f)'{0) = 0. Figure [3] shows both the approximated analytic and numerical (potential included) 
solutions. The approximation becomes worse as the radial distance grows. 
For the parametrized model (PM) : ppM{,f) ^'^ 



= / pPMir)r^dr, 

Mpi 



Mpi 3- a 

where 2-^1 is the hypergeometric function. 



P sH0 2 n , L^ + Ci, (33) 
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Figure 3: The chameleon profile in the ISO DM halo of U5005 galaxy for /3 = 1.69 x 10"''. 
The analytic and numerical chameleon solutions are represented in red line and black line 
respectively. 

Taking limit r — )■ and assuming 0' — )■ oo slower than 1/r^ in this limit, we obtain 
Ci = 0. Then 

= ' , (34) 

Mpi 6 — a 

or 

r/r. 



= 77 o : — • (35) 



13 ri-°r>o(l + ^)°-^2Fi(3 - a, 1, 4 - a, ^) 
Mpi 3- a 

The value of 0'(O) for each a can then be obtained, 



Case I : a < 1 

0'(O) = 0. (36) 

Case II : a = 1 (NFW) 

0'(O) = (37) 

ZIVlpi 

Case III : 1 < a < 2 

(j)'{0) = oo. (38) 
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The corresponding solution of the chameleon scalar field in the PM dark matter halo is 







Mp, 3- a V 2 - a 



-2-Fi(2 — a, 3 — a, 4 — a, —- 



+ 0(0), 



(39) 



or 



(r) = ^^f^(l + ^) 



"-22Fi(2-a,l,4-a, 



r/vs 



1 + r/rg 



+ 0(0). (40) 



where a must be less than 2. Figure H|5|6] show both approximated analytic and numeri- 
cal (potential included) solutions for the PM DM profiles, the differences are hardly visible. 
Figure [7] summarizes the numerical solutions for the PM model. The behaviour of the 
boundary value 0'(O) confirms the analytic results. 

As a summary, we have shown that depending on the DM profile, the chameleon so- 
lution canNOT arbitrarily satisfy the boundary condition 0'(O) = 0. The NFW and the 
parametrized profiles with a > 1 {a must be less than 2 for positivity of the chameleon 
distribution function) cannot satisfy this fiat boundary condition as shown above. 
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Figure 4: The chameleon profile in the PM DM halo of U5005 galaxy for a = 0.2 and 
P = 2.24 X 10~^. The analytic and numerical chameleon solutions are represented in red line 
and black line respectively. 

Now we want to derive a bound on the coupling /3 for the nonsingular solutions which 
require at the galactic edge to match (j)min of the universe. The value of (pmin is determined 
from the observed value of the dark energy through the self-interacting potential of the 
chameleon. Using the analytic formula, we can calculate the maximal matter-chameleon 
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U5005 PM g=0.5 Halo, /3 = 1.96 x 10"'^ 
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Figure 5: The chameleon profile in the PM DM halo of U5005 galaxy for a = 0.5 and 
P = 1.96 X 10"''. The analytic and numerical chameleon solutions are represented in red line 
and black line respectively. 



U5005 PM a=OJ Halo, /3 =1.76 x 10"^ 
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Figure 6: The chameleon profile in the PM DM halo of U5005 galaxy for a = 0.7 and 
/3 = 1.76 X 10"''. The analytic and numerical chameleon solutions are represented in red line 
and black line respectively. 
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Figure 7: The analytic chameleon profiles in the PM DM halo of U5005 galaxy for /3 = 
1.76 X 10~^. The red lines represent the chameleon profile for < a < 1 and the black lines 
represent the chameleon profile for 1 < a < 2. 



coupling of this solution by integrating Eqn. ( ETj) . 



> 



/3 



dr—M{r). 



With 



given by Eqn. ([6]), the maximum Pmax is found to be 



(41) 



1 

n + 2 



47rM 



pi 



r 

Jo 



dr 



M{r) 



71+1 

ri+2 



(42) 



This crucial limit on the matter-chameleon coupling is originated from the non-singular 
boundary condition 0'(O) = 0. As the object gets more and more massive and substantial, 
the maximum value of /3 decreases accordingly. For galaxy U5005 with ISO and NFW DM 
profile, the values of Pmax — 1.69, 1.76 x 10"'^ respectively. 

For the self potential of the form V{(j)) = M'^(l + /^(M/0)"), we can make substitution 
M^"*"" — )■ /iM^"^" to obtain the limit containing two parameters as 



< 

i-NJ 



pi 



Pc 



1 

n + 2 



pi 



r 

Jo 



dr 



M{r) 



n + l 
n + 2 



(43) 



For n = 1, the value of fi is constrained by the LLR experiment and cosmo logical conditions 
to be smaller than 10^ [19j, resulting in the upper limit /3 < W^^^l^max for f^max given in 
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Eqn. fl42p . This is roughly 46 times larger than the original inverse-power-law self potential 
case. For NFW and ISO DM profile of galaxy U5005, the upper limit corresponds to about 
/3 < 10-^ 

Among the three classes of solutions mentioned earlier, only the one satisfying non- 
singular boundary condition was previously considered physically relevant. However, as 
demonstrated above, the non-singular boundary condition canNOT be satisfied for arbitrary 
substantially massive object such as the galactic DM halo. At the time of structure formation 
when the scalar field obtained their vev as the minimum of the effective potential, there is 
no way the chameleon could know which value the matter-chameleon coupling /3 should be. 
This value should have been fixed by some theory at the high scale. For generic situation 
during structure formation, the chameleon should thus develop a profile with singular point 
at the center r = 0. There is no physical reason to prevent this class of solutions for the 
galaxy. 

3.2 Chameleon solutions with singular point 

In general situation, the chameleon could develop singular profile at r = when 7 < Mq. 
For sufficiently large /3, the chameleon solutions with < 7 < Mq yield mostly the same 
results as the case 7 = 0. Without loss of generality, we will therefore consider the profile 
with 7 = and numerically solve for the profile of chameleon within the dark matter halo 
using the boundary condition (pirmax) = 0min, <P'{^max) = for r^ax defined to be the radial 
distance where pmijmax) — Poo, the average density of the universe. We will set n = 1 for 
the chameleonic self-interacting power-law potential V^(0) = and briefiy discuss the 

negligible changes for other n in Section 5. The chameleon profiles for the NFW and ISO 
dark matter halo for some LSB galaxies are shown in Fig. [HI 

The density of chameleon turns out to be only about 10~^ of the DM density and thus 
no distinguishable effects of the chameleon could be seen in the gravitational lensing. The 
gravitational contribution of the chameleon pressure is also negligible. Regardless of its 
smallness as a gravitational entity, the chameleon pressure gradient induces appreciable 
force on the DM halo in the form of the fifth force. 

The chameleon changes from dark energy to matter as the density of the DM increases 
towards the galactic center. This is shown by the ratio w = P^/ p^j, as in Fig. [HI As it changes 
into some form of exotic chameleonic matter with w = 1, the acceleration of the test object 
due to the chameleonic fifth force increases as shown in Fig. [9] 

4 Rotation Curves 

The fifth force induced by the chameleon on the dark matter has variation along the radial 
direction according to Eqn. f [20|) . The circular velocity is then reduced when the gradient of 
the chameleon is negative. The circular velocity which includes the acceleration due to the 
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U4325 , ISO Halo U337 1 , ISO Halo 

5. X 10"^ \ ^ ] f — ■ 




500 1000 1500 2000 200 400 600 800 10001200 
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Figure 8: The chameleon profile and the ratio w — P^/ oi the chameleon in the ISO DM 
halo of U4325, U3371 galaxy for varying ^Q-^-e — 2, 4, 6, 8, 10 (from bottom to top graph). 
The chameleon field and w increase with the coupling (3. 




R(kpc) R(kpc) 

Figure 9: The acceleration from the chameleonic fifth force in the ISO DM halo of U4325, 
U3371 galaxy. The acceleration increases with the coupling /3. 
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fifth force from the chameleon can be written as the following 



GM{r) (5r dcj) 

The accumulated mass of the dark matter halo depends on which mass model we are using. 

The NFW DM profile emerged as the universal density profile for the dark matter halo 
from the simulations of ACDM model. The density profile depends on the critical density 
required to flatten the universe and therefore it is determined uniquely by the Hubble pa- 
rameter H. The profile correctly gives a constant rotation velocity at large radii for most 
galaxies. However, recent observation of the rotation curves from the late-type LSB galaxies 
reveals that the center of these galaxies contains an approximately constant density core 
better parametrized by the ISO DM profile [21J. Among the LSB galaxies in Ref. pT] . 
U4173, U4325, U3371, DD0185, DD047, DD064, U1281, DD052, IC2233 are the galaxies 
which the NFW profile cannot be used to fit with the observational data of the rotation 
curves around the core region without stretching it to unrealistic parameters, namely too 
large asymptotic rotation velocity at large radii. On the other hand, the ISO profile can fit 
well to most LSB galaxies implying that these galaxies might actually have a constant core 
region. 

The contrast between simulation results from the ACDM model and observational data 
is known as the core-cusp problem (see Ref. [22] for an excellent review). While the ACDM 
model supports a cuspy rotation curve, observational data of the late-type LSB galaxies 
prefer the mass model with a constant core. To explore the effects of the chameleonic fifth 
force on each kind of DM halo, we consider two popular mass models, NFW and ISO, and 
finally in a parametrized (PM) model. 

4.1 NFW halo 

For the NFW profile, the DM density is given by 

/'NFw(r) = ^ /° 2 > (45) 

a V a) 

where po and a are the characteristic density and the scale radius of the halo respectively. 
Although we will use the NFW profile and effective acceleration from the chameleonic fifth 
force to calculate the circular velocity, we must know the characteristic density of the NFW 
profile. In order to obtain the po of NFW profile, we consider the accumulated mass 

/ X ..2 / ln(l + (cr/raoo)) - (cr/rapp)/ (1 + {cr/r2m)) \ ^200 
M,,^{r) = [ ln(l + c)-c/(l + c) ) ^' ^^'^ 

leading to the circular velocity in the absence of the chameleon as 



-Mr)=v-.»,,/- '°!\!r'"'f;;f ■ (47) 

ln(l + c) — c/(l + c) 
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where V200 is the circular velocity at virial radius (r2oo) , x = r/r2oo is the radial distance in 
unit of the virial radius and c = r 200 1 a is the concentration parameter. 
The characteristic density of the NFW profile is then given by 



Po 



'^200 



A-kGo? Vln(l + c) -c/(1 + c) 



(48) 



And from the virial mass 
respectively: 



"^200)5 we obtain the virial radius and the characteristic radius 



f 200 



^200 

mi' 
f 200 



(49) 
(50) 



where H is the Hubble parameter. 

For the late-type low surface brightness (LSB) galaxies, dominating gravitational element 
is the dark matter halo. It is thus the most efficient to study effects of chameleonic fifth 
force on the rotation curves of the LSB galaxies. Using values of V200, c from Ref. |2T] and 
H = 72 km/s/Mpc, we can produce the rotation curves for certain LSB galaxies shown in 
Fig.[Tni 



U5005 NFW Halo, 




DD0189 NFW Halo, 



= 1,1.5,2 




Figure 10: Rotation curves of U5005, DD0189 galaxy around the core region for varying 
p. The red lines represent the rotation curve of the galaxy with NFW profile without the 
chameleon. 



4.2 ISO Halo 

For a spherical pseudo-isothermal halo, the DM density profile is assumed to be 



Pisol^^j 



Pc 



2 ' 



(51) 
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U5005 ISO halo, 




DD0189 ISO halo 



= 2,4,6 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 14 
R(kpc) 




1.0 1.5 2.0 
R(kpc) 



Figure 11: Rotation curves of U5005, DD0189, U4325, U3371 galaxy in the core region for 
varying /3. The red lines represent rotation curves of the galaxy with ISO profile without 
the chameleon. 
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where pc and Rc are the central density and core radius of the halo respectively. The rotation 
curve from this density profile in the absence of the chameleon is 



Vc{r) = ^/47rGpei?2 (^l-yarctan(-^)^. (52) 

The best-fit parameters pc,Rc and processed rotation curves of some LSB galaxies are given 
in Ref. [21]. We use these parameters to generate rotation curves with varying chameleon- 
matter coupling P as are shown in Fig. [TTl 



4.3 The parametrized model 

We can generalize the density profile of the NFW model to be 

ppM(r) = — ^33^, (53) 

where a is a parameter which takes value of 1 for the NFW profile. This profile correctly 
reproduces p oc in the large radius limit. The circular rotation velocity in the absence 
of chameleon is then 



—a 



^c(r) = W^^Hj-l 2^,(3 -a, 3 -a, 4- a, -r/r,) (54) 



3 — a \r 

where 2-^1 is the hypergeometric function. Rotation curves of U5005, DD0189 with a 
0.2, 0.7 and U4325, U3371 with a = 0.2 for varying /3 are shown in Fig. | 



4.4 Analytic approximation of the rotation curve 

An excellent approximation of the rotation curve can be obtained analytically according to 
Eqn. (|23|) . Demanding that (f)'{rmax) = {i.e. 7 = 0) at the galactic boundary, we obtain 
{(j)'r'^\r=o) = ~ i-fu ^ -^i^max)- The analytic approximation of the acceleration from the fifth 
force is then given by 

The galactic mass Mq is again defined to be M{rmax)- The force pushes outwardly when 
the chameleon gradient is negative. It effectively reduces the rotation velocity of matter in 
the galaxy. For a generic DM halo which is far from undergoing a gravitational collapse, 
B[r) ~ 1. The acceleration from the fifth force is thus determined mostly by the mass profile 
and it is proportional to Substituting into Eqn. (jUj), the resulting rotation curves agree 
extremely well with the numerical results. 
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U5005 a = 0.7 Halo, 



J_ 
10- 



= 0.5,1,1.5 



DD0189 a = 0.7 Halo, 



= 1,2,3 





1.0 1.5 2.0 
R(kpc) 



Figure 12: Rotation curves of U5005, DD0189, U4325, U3371 galaxy around the core region 
for varying /3. U4325 and U3371 cannot be fit to the PM model with a = 0.7 without making 
the rotation velocity unrealistically large. The red lines represent the rotation curves of the 
galaxy with the PM profile without the chameleon. 
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5 Results and Discussions 



For NFW profile, we choose to explore the effects of chameleon in U5005 and DD0189 galaxy 
since the data can be fit without assuming too large circular rotation velocities. As we can 
see from Fig. [ini the chameleonic fifth force is larger in the core region. The rotation velocity 
is reduced as the fifth force drives the matter in the outward direction. For U5005, turning 
the value of /9 > 6 x 10"'^ gives the rotation curve cuspier, so much that it misses the error 
bar of the innermost data point. For DD0189, /3 > 1.75 x 10~^ results in the bad fit of the 
innermost data point. For statistical analysis, we use U5005 with the number of degree of 
freedom A^ = 11 — 3 = 8 and establish a constraint on the upper bound of matter coupling 
/3 < 1 X 10~^ at 95 % C.L. DD0189 gives less stringent constraint. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.1 0.2 0.3 0.4 



R(kpc) R(kpc) 

Figure 13: Rotation curves of galaxy U5005, U4325 for the chameleonic self potential V{(f)) = 
with n = 0.1 — 4 (left to right) using NFW, ISO profile. The red lines are the 
rotation curves without the chameleon. 

For ISO profile, we present rotation curves for U5005, U4325, U3371, DD0189 galaxies 
with varying /3. The strongest constraint on the coupling /3 is from U4325 where (3 > 3x 10~^ 
results in a large deviation of the rotation curve from the innermost data point. This is 
shown in Fig. [TTl For statistical analysis, we use U4325 with the number of degree of 
freedom A^ = 16 — 3 = 13 and establish a constraint on the upper bound of matter coupling 
P < 1 X 10~^ at 95 % C.L. Other galaxies give less stringent constraints. 

The PM mass model can fit with the rotation curve data better than the NFW in the 
core region for a < 1 since the cusp is less steep. For a > 1, the quality of fit is worse, 
therefore we consider the PM mass model only in a < 1 cases. This is shown in Fig. [T2j 
Again, the chameleonic fifth force makes the cusp steeper, resulting in the worse fitting. For 
U5005, DD0189, large deviations of the rotation curve from the innermost data point occur 
when /3 > 1, 2 X 10~^ for a = 0.7 respectively. For statistical analysis, we use U4325 with 
the number of degree of freedom = 16 — 3 = 13 and establish a constraint on the upper 
bound of matter couphng /3 < 5.4 x 10"^ at 95 % C.L. for a = 0.2. For a = 0.7, U5005 is 
used with A^ = 11 — 3 = 8, the upper bound on the matter coupling becomes /3 < 9 x 10^'^. 
The strongest constraints for each DM halo are summarized in Table [TJ Numerical studies 
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LSB galaxy 


upper bound on f3 at 95 % C.L. 


U5005 (NFW) 


6 X 10-3 


U5005 (ISO) 


2 X 10~3 


U5005 (PM a = 0.2) 


6 X 10-3 


U5005 (PM a = 0.7) 


9 X 10-3 


DD0189 (NFW) 


1.75 X 10-2 


DD0189 (ISO) 


4.8 X 10-3 


DD0189 (PM a = 0.2) 


1.75 X 10-2 


DD0189 (PM a = 0.7) 


1.85 X 10-2 


U4325 (ISO) 


1 X 10-3 


U4325 (PM a = 0.2) 


5.4 X 10-3 


U3371 (ISO) 


2.7 X 10-3 


U3371 (PM a = 0.2) 


9.5 X 10-3 



Table 1: Constraints on the matter-chameleon coupling constant from the LSB galaxies. 

show that the rotation curves change very slightly (cuspier for larger n) for the self potential 
with 72 = 0.1 — 4 as shown in Fig. [131 The constraints are thus mostly the same for the 
power-law potential with n = 0.1 — 4 (only fractionally stronger). 

Observe that the rotation curves are made cuspier for every galaxy in the presence of 
the chameleon, regardless of the mass model of the DM halo. The parameters of each mass 
model need to be fit to the rotation curve in the large radii region and they will be the 
same as the best fit in the absence of the chameleon. Deviation of the rotation curves in the 
presence of the chameleon will appear in the core region where the effect of the chameleonic 
fifth force is the most distinctive. 

6 Conclusions 

Generically, a chameleon scalar field is designed to interact with matter through a conformal 
coupling so that its mass varies with the matter density locally. The chameleon becomes 
massive when the local density is high and the modification of gravity from the modified 
Yukawa potential is thus negligible. The chameleon mechanism makes the scalar field evade 
the constraints on the violation of Equivalence Principle and the fifth force from the Earth 
bound experiments. 

In the situation where matter density varies, the chameleon obtains a profile with spatial 
variation, generating a pressure profile within the matter. The pressure gradient force is 
actually the fifth force induced by the chameleon-matter coupling. We demonstrate that 
for a sufficiently large massive object such as a galaxy, the non-singular boundary condition 
leads to stringent but yet unphysical constraint on the matter-chameleon coupling (roughly 
of order /3 < 10"^). We argue that the physical chameleon profile in general inevitably 
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develop a singularity at the center. The chameleon profile with central singularity leads to 
a significantly large chameleonic fifth force close to the center region. In a presence of the 
galactic DM halo, the fifth force becomes appreciably large so that in the core region of the 
galaxy, the force significantly reduces the circular rotation velocity of the galaxy, resulting 
in a cuspier rotation curve. 

We investigate the effects of the chameleonic fifth force on the rotation curves of certain 
late-type LSB galaxies using the NFW, ISO, and parametrized mass profiles. Constraints on 
the upper bound at 95 % C.L. on the chameleon-matter coupling are established in Tabled] 
The constraints could be as stringent as /3 < 1 x 10~^ for the ISO mass model. For a 
parametrized mass model with a = 0.7, 0.2, the constraints are not as strong, /3 < 5— 9x 10~^. 

Analytic approximation of the rotation curve in the presence of the chameleon is derived 
with great accuracy. The chameleonic fifth force is determined mostly by the mass profile of 
the DM halo and it is proportional to the matter coupling The force is directed outward 
since the chameleon gradient is negative. The force is strongest in the core region and the 
rotation curve is altered the most in the region. The chameleon therefore generates cuspier 
rotation curves for every mass model of the DM halo. 
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A Derivation of the chameleonic fifth force 

From the chameleon-matter action 



after varying with respect to the field, 0, we obtain an equation of motion of the chameleon 

as 



The above equation must be consistent with the equation from V^T^y = 0. Since the 
chameleon scalar field couple with matter, then 



S 




v20 = v^,^-a^r;(™). 



(56) 




(57) 



The energy-momentum tensor of the scalar field is 
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Thus 



where 



Then 



9^0(v^5,0) - 9>v.9,0 = 9^0(a'^a,0-^7'^'^r^,9A0)- 9X9,9,0 -r^,9A0), 

0. 



-9°0r^,9A0 + 9-0r;i,9A0, 



Therefore 



VmtW = (V20-V^,0)9,0. 
From the equation of motion of the chameleon (Eqn. (!56l) ). we obtain 

= a^pdycj). 

From Eqn. (IFTIl . the covariant derivative of the matter is then 

V^T^::) = -«^p9.0. 
For non-relativistic fluid and approximately flat spacetime, the LHS becomes 

We approximate the matter to be pressureless while the gradient has only the radial direction 
due to spherical symmetry ■ V = 0). Therefore 

p(9if*) = -a^pdi(j), 

then 

a = — a(^V0. 

This is the acceleration due to the fifth force from the chameleon acting on the matter. 
Therefore, if we consider only the effects on the matter in presence of gravity , we will obtain 
the circular velocity of matter within the galaxy as 



GM{r) 

'^c{r) = \ ha^(9r0)r. 
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